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Summary 


The objective of this effort is to develop an analytical model for the coupling of 
active noise control (ANC) piston-type actuators that are mounted flush to the inner and 
outer walls of an annular duct to the modes in the duct generated by the actuator motion. 
The analysis will be used to couple the ANC actuators to the modal analysis propagation 
computer program for the annular duct, to predict the effects of active suppression of fan- 
generated engine noise sources. This combined program will then be available to assist in 
the design or evaluation of ANC systems in fan engine annular exhaust ducts. 

An analysis has been developed to predict the modes generated in an annular duct 
due to the coupling of flush-mounted ring actuators on the inner and outer walls of the 
duct. The analysis has been combined with a previous analysis for the coupling of modes 
to a cylindrical duct in a FORTRAN computer program to perform the computations. 
The method includes the effects of uniform mean flow in the duct. The program can be 
used for design or evaluation purposes for active noise control hardware for turbofan 
engines. 

Predictions for some sample cases modeled after the geometry of the NASA Lewis 
ANC Fan indicate very efficient coupling in both the inlet and exhaust ducts for the m = 6 
sp inning mode at frequencies where only a single radial mode is cut-on. Radial mode 
content in higher order cut-off modes at the source plane and the required actuator 
displacement amplitude to achieve 1 10 dB SPL levels in the desired mode were predicted. 
Equivalent cases with and without flow were examined for the cylindrical and annular 
geometry, and little difference was found for a duct flow Mach number of 0. 1 . 

The actuator ring coupling program will be adapted as a subroutine to the 
cylindrical duct modal analysis and the exhaust duct modal analysis. This will allow the 
fan source to be defined in terms of characteristic modes at the fan source plane and 
predict the propagation to the arbitrarily-located ANC source plane. The actuator 
velocities can then be determined to generate the anti-phase mode. The resulting 
combined fan source/ANC pressure can then be calculated at any desired wall sensor 
position. The actuator velocities can be determined manually or using a simulation of a 
control system feedback loop. This will provide a very useful ANC system design and 
evaluation tool. 
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1. Introduction 

1.1 Program Objective 

The objective of this effort is to develop an analytical model for the coupling of 
active noise control (ANC) piston-type actuators that are mounted flush to the inner and 
outer walls of an annular duct to the modes in the duct generated by the actuator motion. 
The analysis will be used to couple the ANC actuators to the modal analysis propagation 
computer program for the annular duct, to predict the effects of active suppression of fan- 
generated engine noise sources. This combined program will then be available to assist in 
the design or evaluation of ANC systems in fan engine annular exhaust ducts. 


1.2 Previous Effort 

In a prior analysis 1 , a model was developed for the coupling of ANC sources to a 
cylindrical duct, more representative of a turbofan inlet. The reader is referred to that 
paper for background on prior research and motivation for the application of ANC 
systems to aircraft engines. 

The previous analysis did not account for the effects of mean flow, which is 
included in this analysis. The present study is based in large part on the development 
presented in the previous paper, extending the analysis to include annular duct geometry 
with ANC actuators on both inner and outer walls and the effects of uniform mean flow. 
The results apply to both cylindrical and annular ducts, with or without mean flow. The 
extension to annular ducts provides a first approximation to simulating a fan exhaust 
system. 


1.3 Duct Geometry 

Figure (1). shows the duct geometry and the mounting positions for the actuator 
rings. The radial variable is r, the angular variable is 0, and the axial variable is z. The 
duct inner and outer radii are given by n and r 2 , respectively. The actuator ring has axial 
extent L, and the zero of the axial coordinate system is taken as the upstream end of the 
actuator ring The mean flow velocity is denoted as U, and the flow is positive when 
moving in the direction of the positive z-axis (exhaust duct flow will be positive, inlet duct 
flow will be negative). 
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Figure (1). Annular geometry for ANC actuator mounting. 

1.4 Assumptions 

The analysis is subject to a number of idealizations and simplifying assumptions. 
Among these are: 

1 . The duct geometry is axisymmetric and is assumed to be constant in inner and outer 
radii along the full length of the duct. 

2. The mean flow is assumed to have a uniform profile, that is, contains no shear layers. 

3. The actuators are rectangular in shape, and behave like rigid pistons, achieving a 
volume velocity equivalent to a piston mounted flush to the wall. 

4. It is assumed that there is no interaction between the propagating fan-source- 
generated mode passing the actuator rings and the pistons, that is, the pistons are 
assumed to have infinite internal impedance. 

The program is meant primarily for computing the case in which only one radial 
mode is cut-on. The actuator velocity that generates the first radial mode that matches 
the fan source mode is found. The complex ratio of inner actuator velocity to outer 
actuator velocity is given, and this ratio is maintained throughout the computation as the 
outer wall actuator velocity is scaled. The program solution determines the pressure 
(modal) spillover into higher-order cut-off modes at the actuator source plane. Although 
these modes do not propagate, the smaller their value the more efficient is the coupling of 
the actuators to the desired mode. 

The program can also be used for the case with multiple cut-on radial modes, but 
the computation attempts to match only the first radial mode to the input SPL level. The 
given inner actuator velocity to outer actuator velocity ratio is maintained throughout the 
computation. The mode coefficient obtained for the second radial mode will depend on 
the velocity chosen for the inner actuator, but the program does not attempt to find the 
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inner actuator velocity that best matches both fan source mode coefficients 
simultaneously. This, in fact, may be a very difficult problem to solve. 

1.5 Approach 

The approach will be to use a classical modal solution to the wave equation in an 
annular duct with non-homogeneous boundary conditions. A transformation of variables 
method is applied that transforms the homogeneous differential equation with non- 
homogeneous boundary conditions to a non-homogeneous differential equation with 
homogeneous boundary conditions, which is more tractable. A modal expansion in 
hardwall duct modes is used to solve the radial variation and the axial variation is solved 
using infinite Fourier Transforms. This is similar to the approach in Reference 1, with the 
extension to more complex annular geometry and the inclusion of the effects of mean 
flow. 


1.6 Summary of Results 

An analysis has been developed that predicts the modes generated in an annular 
duct due to the coupling of flush-mounted ring actuators on the inner and outer walls of 
the duct. The analysis has been combined with a previous analysis for the coupling of 
modes to a cylindrical duct in a FORTRAN computer program which is operational and 
partially validated. The method includes the effects of uniform mean flow in the duct. 
The program can be used for design or evaluation purposes for active noise control 
hardware for turbofan engines. 

Predictions for some sample cases modeled using the geometry of the NASA 
Lewis ANC Fan indicate very efficient coupling in both the inlet and exhaust ducts for the 
m = 6 spinning mode at frequencies where only a single radial mode is cut-on. Radial 
mode content in higher order cut-off modes at the source plane and the required actuator 
displacement amplitude to achieve 1 10 dB SPL levels in the desired mode were predicted. 
Equivalent cases with and without flow were examined for the cylindrical and annular 
geometry, and little difference was found for a duct flow Mach number of 0. 1 . The use of 
the inner ring of actuators to improve coupling efficiency to a single radial mode and to 
couple to two radial modes simultaneously was examined. 
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2. Analysis of ANC Actuator Coupling to an Annular Duct 

2.1 Wave Equation for the Acoustic Pressure 

The partial differential equation for the acoustic pressure in a duct with uniform 
mean flow is given by 




— + U — 
dt dz 


2 

I P 


= c 2 V 2 p 


0 ) 


where 


t 

z 

P 

c 

U 


time 

axial variable 
acoustic pressure 
speed of sound 

mean flow velocity in the z- direction 


Using the Laplacian operator in a cylindrical coordinate system and postulating harmonic 
time dependence of the form 


P(t) = e_i<Bt 


(2) 


where 

i = V-T 

co = circular frequency, = 27rf 

f = frequency, Hz., 


we can write Equation ( 1 ) in the form 

= \ 2 


1 


-ico + U — 

dz 


1 8 
P = “T- 

r dr 
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\ 1 d 2 p c^p 

+ — s- — £ + — 


or 


r 2 50 2 


dz 


(3) 


where 

r = radial variable 

0 = circumferential variable. 

The periodic boundary condition in 0 allows us to separate the 0-dependence of 
the solution in the form 
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P(r,z,0) = p(r, z)e : 


(4) 


im0 


where 


m = spinning mode order, any positive or negative integer or zero 

When Equation (4) is substituted into the partial differential equation, this further reduces 
the form to 


1 8 _ 
r dr 


r 



m 2 w dp 

-^-p + 2ikM-^ + 

r z OZ 


(l + M 2 ) 


^f + kVo 

dz 2 


where 

k = wavenumber, = co/c (constant) 

M = mean flow Mach number, = U/c (constant) 


(5) 


This is the homogeneous partial differential equation in the radial and axial 
variables. The solution for p in r and z must satisfy this equation and the boundary 
conditions applied by the actuators at their locations on the inner and outer walls of the 
duct. 


2.2 Boundary Condition with Slip Flow 

The actuators will be modeled as rigid pistons of rectangular shape that vibrate 
with infinite internal impedance at frequency co. On the actuator surface, the wall has a 
radial velocity component denoted by v w . The boundary condition, which must include the 
effects of slip flow at the wall surface, will be written in terms of the requirement that the 
normal component of acoustic particle displacement be continuous at the wall (it should 
be noted that the physically incorrect assumption of slip flow leads to the anomaly of 
different boundary conditions depending on whether particle displacement or particle 
velocity continuity is assumed 2 ). 

The normal velocity of the actuator surface at the wall is related to the surface 
displacement for harmonic time dependence as 

v w = -ioC ( 6 ) 

where C, denotes either the actuator surface displacement or the acoustic particle 
displacement at the wall. The acoustic particle velocity at the wall is written in terms of 
the acoustic particle displacement using the convected derivative, 
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where v r is the radial component of the acoustic velocity. 

The radial component of the acoustic momentum equation is given by 


— = ipcov r - pU 
dr 


dz 


( 8 ) 


where p is the ambient air density. The assumption of continuity of particle displacement 
at the wall requires that the actuator displacement be the same as the normal component 
of acoustic particle displacement. Combining Equations (6) and (7) gives the relation 
between actuator normal velocity and the acoustic particle velocity normal to the wall. 


v r = 


1 + i 


.M d 


k dz 


(9) 


Substituting the expression for the acoustic velocity at the wall into the momentum 
equation gives the boundary condition for the pressure at the wall 


op 

dr 




( 10 ) 


The actuator velocity at the wall can be Fourier-analyzed into its various spinning 
mode components in the circumferential direction, which will depend on the functional 
form of the actuator velocity in 0 


v,(6)= Xv„«*' (11) 

m = -oo 

where Wxm is the m-order component of the actuator circumferential velocity profile at the 
wall. This expansion and a method to determine the v^ to generate a preferred m-order 
were the considered in Kraft and Kontos\ For this analysis, we assume that only a single 
m-order is under consideration, and that the m-order component of the actuator wall 
velocity variation has been determined and is given. 

The actuator velocity at the inner and outer walls can be expressed in terms of the 
m-order spinning mode component and the axial velocity profile as 


v l» = v lm h (z) 


( 12 ) 



and 


v 2w = v 2m h (z) 


(13) 


where subscript 1 refers to the inner wall, subscript 2 refers to the outer wall, and h(z) is 
the dimensionless axial functional dependence of the actuator velocity. 

The boundary condition equation at the inner wall can then be written 


dp 

dx 


iM d_ 
k dz 


A lm h(z) 


atr = r ls 0<z<L 


(14) 


and the boundary condition at the outer wall can be written 

A 2m h(z) atr = r 2 , 0<z<L 


ap 

dr 


f iM d ^ 2 

1 + 

k dz 


(15) 


where we define the m* order coefficient of the wall velocity at the inner wall as 

Ai m = -ipckvjm 

and the m 01 order coefficient of the wall velocity at the outer wall as 

A 2m = -lpckv 2m 


(16) 


(17) 


We have made an implicit assumption here that leading and trailing edges of the inner and 
outer rings align at the same axial stations. Later it will be shown that this is necessary to 
simplify the analysis. Rings that do not align axially can be considered separately and the 
generated pressure fields will add linearly. 

We shall also refer to the effective volume velocity of the actuator rings, which is 
defined at the inner wall as 


v im =e me Jv lm h(z)dz (18) 

o 

and at the outer wall as 

V 2m = e ime Jv 2m h(z)dz (19) 

o 
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For the purposes of this analysis, we assume that the duct walls adjoining the 
actuator rings at both the inner and outer walls are hard, so that the boundary condition 
away from the actuator region is given by 


^ = 0 

dr 


at r = q and r = r 2 for z < 0 and z > L 


( 20 ) 


2.3 Transformation to a Non-Homogeneous Differential Equation with 
Homogeneous Boundary Conditions 

The objective is to transform the homogeneous differential equation with non- 
homogeneous boundary conditions into a non-homogeneous differential equation with 
homogeneous boundary conditions. This is accomplished using a three-dimensional 
transformation of variables for the acoustic pressure of the form 

f (r, 9, z) = p( r, 0, z) - g(r, 9, z) (21) 

where g(r,0,z) is chosen to be an otherwise arbitrary function that satisfies the non- 
homogeneous boundary conditions. 

The choice we shall make for g(r,0,z) is a quadratic form in r: 

g(r,0, z) = [B ir + B 2 r 2 ]e im0 h(z) (22) 

where Bi and B 2 are undetermined constants. Then, at r = n, 

^ = [B 1 +2B 2 r,]e ime h(z) (23) 

and at r = r 2 , 

~ = [Bj + 2B 2 r 2 ]e ime h(z) (24) 

We want to show that f(r,0,z) satisfies homogeneous boundary conditions at ri and 
r 2 . The radial derivative of f at r = ri can be written as 

T = ?-7T = f 1 + i Tr] A 2m h(z)-[B,+2B 2 r,]h(z) (25) 

dr or or ( k dzj 
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where Equation (15) for the boundary condition in p and Equation (24) for the boundary 
condition in g have been substituted and the common e un9 factor has been divided out. 
Similarly, the radial derivative of f at r = r 2 can be written as 


<9f _ dp dg _ 

dr dr dr 


. m d 

1 + 1 

k dz 


\ 2 


A 2m h ( z )-[ B i +2B 2 r 2 ]h(z) 


(26) 


To get homogeneous boundary conditions in f, we set 


5f n 

— = 0 at r = n and r = r 
dr 1 


(27) 


Carrying out the derivatives of h(z) from the square of the converted derivative in 
Equations (25) and (26), we can define a function 


tt / n , M h' 
Hm( z ) - 1 + 2i— — 
k h 


mY— 

k J h 


(28) 


where the primes denote differentiation with respect to z. Equations (25) and (26), then, 
are two equations in two unknowns for Bi and B 2 , which can be solved to give 


Bj - H M A 2m - 


Hm(A2hi - A lm ) 


r, -r. 


(29) 


and 


_ ^M(A 2 m A lm ) 

2 ( r 2- r l) 


(30) 


Note that Bi and B 2 are functions of z. With the coefficients Bi and B 2 , the function g is 
completely determined, and the function f satisfies homogeneous (zero derivative) 
boundary conditions. 

It remains to derive the differential equation that f must satisfy. This is found by 
substituting Equation (21) for p into the differential equation. After a fair amount of 
algebra, the result can be shown to be 
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ld_ 

r ar 


af') 

dr J 


m 


■f + 2ikM 


df 

dz 



+ k 2 f = 




A 


Bj +(4 + m 2 -k 2 r 2 )B 2 


J 


h - 2ikMr[B ] + rB 2 ]h' -(l - M 2 )r[B, + rB 2 ]h" 

(31) 


This is the non-homogeneous partial differential equation that must be satisfied by the 
function f, which, by virtue of the particular choice of g, has homogeneous boundary 
conditions. 


For convenience, we shall write the diflferential operator on the left-hand-side of 
Equation (31) as 


i() 


\d_ 

r dr 


6 r 


-^ + 2ikM — + (l + M 2 )-^=- + k 

r - 2 dz \ > A -, 2 


dz 1 


(32) 


The right-hand-side of the diflferential equation can also be written in a somewhat easier to 
grasp form as 


Q(r, z) = Dj - + D 2 + D 3 r + D 4 r 2 
r 


(33) 


where 


0 

N 

III 

IT 

1 

'jiEf 

It 

(34) 

D 2 (z)s(4 + m 2 )B 2 h 

(35) 

D 3 (z) = - 

k 2 h + 2ikMh' + (l - M 2 )h"]B, 

(36) 

D 4 (z) s - 

k 2 h + 2ikMh' + (l - M 2 )h"]B 2 

(37) 


Using Equations (32) and (33), the differential equation can now be written simply 
as 

2(f) = Q(r,z) (38) 
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Note that the “source” term on the right-hand-side of the differential equation 
depends on frequency, Mach number, radius, and the z-dependence of the actuator 
velocity profile. 


2.4 Derivation of the Non-Homogeneous Differential Equation in z 

The solution to the non-homogeneous differential equation is obtained as an 
expansion in a series of eigenfunctions of the radial modes in an annular duct such that the 
expansion coefficients are functions of z that depend on mode order (m,n) 


f( r >z)= X^nn( z ) ( Pmn( a mn r ) ( 39 ) 

n=l 

where n is the radial mode index ordered such that the lowest order radial mode (the first) 
has n = 1, and the annular duct eigenfunctions are given by 

?mn (^mnO — 3 m(^irm0 ^mn^m (^mnO (40) 

where J m is the rn 11 order Bessel function of the first kind, and Y m is the m* order Bessel 
function of the second kind. The eigenvalues a™ are the roots of the boundary condition 
equation for a hardwall annular duct, 

g(p mn( a m rT) =Q a t r = r, and r = r 2 (41) 

dr 

These roots can be found by various means, and are assumed to be known. Root-finding 
subroutines are provided as part of the program that implements this analysis for both 
cylindrical and annular ducts. 

If we define 


Y mi) — ®mn*2 

then the mode coefficient C™ can be written as 


rc* 3 m(y mn) Y mn 3 m+l(Y mn) 

mY m (y 

mn) Y mn^m+l(Y mn) 


(42) 


(43) 


Note that if r } = 0, them C™ = 0. 
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The first step in finding the unknown coefficients fmn is to multiply through the 
differential equation (38) by cp^r and integrate from ri to r 2 . 


J cPmj^f )rdr = J Q(r, z>p m j(a mj r)rdr 


(44) 


Substituting Equation (39) for fin the differential operator, we have 


40 = Z4 f nm( z ) ( Pmn(anm r )] 

n=l 


(45) 


Applying the differential operator to the individual terms in the sum, we get 

f \ _2 


4 f mn<Pmn] = f rr 


l_e_ 

r dr 


V 


dcPmn ] _ m ^mn , i^2 

*3 2 ^ K ^ mn 

dr j r z 


+ 2ikM^<p ran + 

dz 


(l-M 2 ) 


oz 2 


■<Pr 


(46) 


The first two terms of the term in square brackets on the right-hand-side of Equation (46) 
can be recognized as the Bessel eigenfunction equation for the radial mode in a cylindrical 
or annular duct, so that we can (without proof) state the following relationship: 


\_d_ 
r dr 


f 


<9<P r 


■ mn |_ g 1 9_mn _ 2 

11 ? u mnYr 


V 


dr 


(47) 


Making these substitutions. Equation (44) becomes 


1 12 

(l-M 2 )f"j + 2ikMf^j +(k 2 -a 2 mj )f mj |Q(r,z)(p mj (a mj r)rdr (48) 

r, 


where 


"Z 

N mj = J <Pmj<Pm/ dr 


(49) 
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is the mode normalization integral, where we have made use of the orthogonality 
properties of the eigenfunctions, 

r 2 

J <Pnrn<Pmj rdr = 0 for n *J- ( 5 °) 

r i 

and where Q(r,z) is given by Equations (33) through (37). 

The integration over r on the right-hand-side of Equation (48) breaks into four 
separate integrals, and can be written as 


J Q(r, z)cp irm (a mn r)rdr = D x (z&P + D 2 Hr + D 3 Hr + D 4 H^ m ( 5 1 ) 

r l 

where 


*2 

Hr = j <Pmn( a mnr) dr 

r i 

*2 

HP = J ( Pmn( a mn r ) rdr 

r l 

H r = J 9mn( a nm r ) r2dr 

r l 

and 

^2 

Hr = J <Pmn(anmr) r3dr 

ri 


(52) 


(53) 


(54) 


(55) 


These integrals must be calculated numerically, but will be constants depending only on 
mode order when computed. The integration subroutines are described briefly in 
Appendix A. 

We now transform the differential equation in z to a form that explicitly shows the 
dependence on the actuator velocity function h(z) on the right-hand-side. The result is 
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( 56 ) 


(l - M 2 )f " + 2ikMf^ + (k 2 - a 2 mj )f mj = Qffi, h + Qg>h' + Qffih' 


(h') 2 


(5) n (6) (h") 2 


0 (4)i£f +0 (5)ii^_ + 0 
Vmn , T Vmn T Vmn 


where 


Qmn = (m 2 - l)BAjHo n + (4 + m 2 )BA 2 H 1 mn - k 2 BA 1 H^ ,n - k 2 BA 2 Hf n (57) 


Qmn = (m 2 - l) 2i — BA]!!™ +(4 + m 2 )2i— BA^” 111 -^MkBAjHf 1 - 


M 


4iMkBA 2 H 


mn 


(58) 


0^ = 
Vmn 


(m 2 -l)^yl BA ] Hr-(4 + m 2 )^Yj BAjHj™ - (l - 2M 2 )BA ] Hf 1 


(l-2M 2 )BA 2 H3 in 

(59) 

(60) 


Qffi=4M 2 [BA 1 Hr-BA 2 H 3 mn ] 

(61) 

Q ( m ^ =-^-(l-2M 2 )[BA 1 Hr - BA 2 H“ n j 

(62) 

M 2 (l - M 2 ) , , 

QSS = —^73 -[bA,H5“ - BAjHf] 

K 

(63) 

and where 


BA, = A 2m 2m lm r 2 

r 2 _r l 

(64) 
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( 65 ) 


__ ^2m ^lm 
2 2 ( r 2- r l) 


The coefficients are constants, so that all the z-variation on the right-hand-side of 
Equation (56) is contained in the h-function and its derivatives. This completes the 
derivation of the ordinary differential non-homogeneous differential equation in z, which 
will now be solved. 


2.5 Fourier Transform for z-Dependence 

We solve Equation (56) using the method of infinite Fourier Transforms. Since we 
have not specified a particular form for the actuator velocity profile function h(z), we will 
only be able to express its Fourier transform symbolically at this point. The last three 
terms on the right-hand- side of Equation (56), which contain products and ratios of the h- 
function and its derivatives, are particularly troublesome. 

The Fourier Transform of the f-function will be written as 


F(K)s?[f(z)]= Jf(z)e +iKZ dz (66) 


The relationships for the Fourier Transform of derivatives of the f-function are 

•?[f'(z)] = -ixF(K) (67) 


and 


?[f"(z)] = -K 2 F(K) 
and the inverse transform is defined as 


( 68 ) 


f(z) = J e 1K2 F(K)dK 
271 J 


(69) 


Note that this form of the Fourier Transform pair is the complex conjugate of that found 
in many reference works. The choice of sign on the exponent in the integrand is arbitrary. 

We will denote the Fourier Transform of the h-function as 

H( k) = <7[h(z)] 
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Taking the Fourier Transform of Equation (56), using the relationships for Fourier 
Transforms of derivatives, and solving for the Fourier Transform of f™, we get 


F mn( K ) = 


fo (1) -kO (2) -k 2 0 (3) ' 

xmn ^Vmn K xmn 

%)+qSX+Q™h 5 +q<2h 6 


(k 2 “ a mn) 

+ 21cMk - 1 

[l-M 2 ] 

|K 2 ' 



(70) 


where FL», H 5 , and H6 are the Fourier Transforms of the functions involving ratios of h(z), 
defined as 




< h ' 2 


(71) 


H< = <? 


h'h " 


(72) 


and 






v h J 


(73) 


There are six inversion integrals to be evaluated: 




Qmn^'( K )e 


— 1KZ 


271 i (k 2 - aL) + 2kMK - (l - M 2 )k : 


■die 


(74) 


G 2 W = -!-j 


-i K Qmn^( K ) e_11 


1KZ 


271 i (k 2 - aL) + 2kMK - (l - M 2 )k : 


die 


(75) 




-K 2 Q ( ^(K)e' 


-1KZ 

dK 


271 i (k 2 - a^i) + 2kMx - (l - M 2 )k 


(76) 


G f ,x 1 J QffiH 4 (K )e- iKZ 

^ J /i 2 ^,2 \ . /i 


2tt ^ (k 2 - «y + 21 cMk - (l - M 2 )k : 


-dx 


(77) 
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(78) 



QmnHg(K)e~ iKZ 

(k 2 - a 2 mn ) + 2kMK - (l - M 2 )k 2 


dK 


and 



QSH.OOe-^ 

(k 2 - a^„) + 2kMK - (] - M 2 )k 2 


dK 


(79) 


The first three integrals will be evaluated using the residue theorem, and the last three 
integrals will be shown to take simple forms for the particular form of h(z) chosen. 

To apply the residue theorem, we need to find the roots in the denominator of the 
integrand, which is the same in all the integrals. We assume that the Fourier transform of 
h(z) possesses no singularities so that all the roots are obtained from the denominator. 


-(i-M 2 )k 2 +2kMK+(k 2 -a^) = 0 

There are two values of k that satisfy this relationship: 

Ji) -kM±^k 2 -(l-M 2 )a 2 „ n 
K 1 - M 2 


(80) 


(81) 


These can be identified as the axial propagation constants for forward and backward 
traveling waves in a duct with uniform flow. In terms of the two roots, we can write the 
denominator of the integrand as 

(l - M 2 )k 2 - 2kMx - (k 2 - i) = -(l - M 2 )(k - k (1) )(k - k (2) ) (82) 


Substituting this into Equation (74) for Gi(z) gives 


G,(z) 


Qj °r %(K)e IKZ dz 
2*(i-M 2 )_ J .(k-kW)(k-k< 2 !) 


(83) 


The inverse Fourier Transform can be evaluated in this case using the residue 
theorem, and the method is considered in detail in Appendix B. The solution for the 
forward-traveling (m,n) mode can be shown to be 
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, Q w 

2 A 2 -(l-M 2 )a; 


where 


-kM+Jk 2 -(l-M 2 )c 


1 - M 

is the axial propagation constant for the forward-propagating (m,n) mode. 
Similarly, the solutions for G2 and G3 can be obtained as 


G 2 (z) = - K mnQ 


+ n (2) 


e ,KmnZ 


fk 2 -(l-M 2 )a; 


G 3 (z) = i(K^ n ) 2 Q ( r 


Q// + 

\ iK 

l_z 

^~ K nmJ e 



-M 2 

)«mn 


We can evaluate the G 4 , G5, and G6 functions if we make some restrictions on the 
form of the ratios of h-functions and derivatives in H4, H5, and H6. If we require that H4, 
H 5 , and H6 have no poles to contribute to the residues within the contour of integration, 
the solution for the G 4 , G5, and G 6 functions can be written in terms of H4, H 5 , and H« as 


G 4 (z) = r 22 " 


iQmn H4 ( Kmn ) g “ 

2 Jk 2 -(l-M 2 )a? 


iQ^ 5 ) 

lv <mn 

h 5 (- 

K mn ) c 

2 

A 2 - 

(l-M 2 )a 2 on 
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(90) 


G 6 (z) = - 


iQmn H6 ( Kmn ) 6 ^ 
2 ^k 2 -(l-M 2 )a 2 


To derive these inverse Fourier transformations, we have had to apply restrictions 
on the form of the actuator velocity z-dependence. We require that the Fourier transform 
of h(z) have no poles in the lower half plane (see Appendix B), and that the Fourier 
transforms of the derivative ratio functions also have no poles in the lower half plane. 
Without these restrictions, it is necessary to return to the original Fourier Transform 
integrals, which may be very difficult to evaluate. 


2.6 Solution for the Acoustic Pressure Field Generated by the Actuators 

The (m,n) mode contribution for the f-fimction that was used to transform the 
differential equation and boundary conditions has been determined and can be written 

f„ (z) = r' [F m (K )] = G, + G 2 + G 3 + G„ + G 5 + G 6 (91) 


writing out the G-terms, this becomes 


fmn( z ) = 


*Qmn K + q(2) + j/ K + ) 2 q( 3 ) 
- ^mnVmn TA ^ K mnj Vmn 


1 %\ 

(~ K ») 

je iK ™ z 

Jv 

Ik 2 -1 



(92) 


1 

2 L 


QL 4 lH 4 (-K^) + Q<lH 5 (-^ n ) + QaH 5 (-Kl)l 6 - — 

yjk -(l-M“)a mrl 


Then, combining Equations (21), (22), and (39), we can express the pressure field 
in the m* order spinning mode generated by the wall-mounted inner and outer actuator 
ring source as 


p(r,e, z) = ||;[f nm (z)(p nm (a mn r)] + [b^z^ + 2B 2 (z)r 2 ]h(z)j 


^im9 


(93) 


The only unknown left to make this a computable formula is the Fourier Transform of the 
actuator velocity profile in z. This, of course, depends upon the choice of function used to 
represent the piston vibration, and will be solved for the particular case of a box function 
in the next section. 
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3. Solution for Box Function Actuator 


3.1 Choice of Box Function Actuator Velocity Profile 

Several choices present themselves for the shape of the axial velocity profile for 
the actuator motion, h(z). The preferred choice will have a known and acceptable Fourier 
transform and will furthermore have a simple Fourier Transform for the three ratio 
functions involving h(z) and its first and second derivatives. It is assumed that the 
circumferential velocity profile has already been accounted for in the spinning mode 
decomposition. 

Three more obvious possibilities are: 

1 . Delta function source - zero axial length but finite volume velocity 

2. Box function source - rigid piston motion with a displacement and velocity 
discontinuity at the edges 

3. Half-sine-wave source - good representation of the lowest simply-supported plate 
vibration mode, displacement and velocity go to zero at the edges. 


The Fourier Transform of the delta function is well known, but some difficulty may 
be encountered in finding the transforms of combinations like 


f 

\ 


g(z> 

5(z) 


2 


5(z) 


and 


S(z) J 


Also, a function of finite extent may be a better physical representation of the actuator for 
the shorter modal axial wavelengths. For these reasons, the delta function was not 
considered at this time. 

The box function, as illustrated in Figure (2) has a constant velocity amplitude 
between 0 < z < L, and represents the motion of a rigid piston actuator. The Fourier 
Transform of this function is given by 


H~ k) 


1 e ~ ilcL -1 
■JT/k -ix 


(94) 


Since the first and second derivatives of h(z) are zero (neglecting edge effects), we have 

H 4 = H 5 = H 6 = 0 (95) 

This is the velocity function used for further illustration by a sample calculation. 
The box function could be an approximate representation of a source of any other shape 
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by choosing the velocity of the piston based on an effective volume velocity for the more 
complex actuator. 


^ Actuator Piston 


Wall 

i 

i 

i 

i 

i 

Actuator 1 ... „ 

Motion I Wa " 

2 = 0 

> 

_l 

|| 

N 




Figure (2) Diagram of box function axial velocity profile for actuator. 

The very realistic half-sine-wave velocity profile has an easily calculated Fourier 
transform, but the transform of the ratio functions may take some additional effort. 
Examination of this case is left for the future, and is recommended as a follow-on task. 


3.2 Sample Computations for Box Function Case 

The computer program FDANC has been written to compute the acoustic pressure 
coupling of wall-mounted rings of ANC actuators in cylindrical or annular ducts. The 
program was used to compute the modes generated in the inlet and exhaust ducts of an 
analytical model that simulates the NASA ANC Fan that is mounted in the NASA Lewis 
AAPL Test Facility. The parameters that are held constant throughout the computations 
are the following: 

Inlet and Exhaust outer radius = 2.0 feet 

Exhaust inner radius = 0.6145 feet 

Flow temperature = 70 degF 

Spinning mode order, m = 6 

Number of radial modes = 6 

Outer wall SPL = 1 1 0 dB 

Axial length of ANC actuator = 4.5 inches 

For this study, both cylindrical and annular ducts were considered, using the above 
parameters Two values of mean flow Mach number were used, Mach 0.0 and 0.1. The 
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annular duct case was run with actuators on the outer wall only, inner wall only, and both 
walls. For the cases with actuators on both the inner an outer walls, the inner wall 
actuator velocity was set to a value based on the ratio of the first radial mode shape at the 
inner wall to that at the outer wall. Two cases with a higher radius ratio (R 1 /R 2 = 0.75) 
were compared. Finally, a case at a higher frequency with two cut-on radial modes is 
presented. A brief outline of the differences among the cases is presented in Table (1): 


Table (1) Differences among sample cases used for computations. All cases have 
same outer radius (R 2 = 2.0) and same spinning mode order (m = 6). 


Case 

Number 

Inner 
Radius 
Ri, ft 

Mach 

Number 

Frequency 

Hz. 

Inner 

Actuator 

Active? 

Outer 

Actuator 

Active? 

1 

0.0 

0.0 

935 

n/a 

Yes 

2 

0.6146 

0.0 

935 

No 

Yes 

3 

0.0 

0.1 

935 

n/a 

Yes 

4 

0.6146 


935 

No 

Yes 

5 

0.6146 

0.1 

935 

Yes 

No 

6 

0.6146 

0.1 

935 

Yes 

Yes 

7 

1.5 

0.1 

935 

No 

Yes 

8 

1.5 

0.1 

935 

Yes 

Yes 

9 

1.5 

0.1 

1350 

Yes 

Yes 


3.2.1 Computer Output for Case 1, Cylindrical Duct, Mach 0.0, m = 6 


This is a basic run for the inlet duct to illustrate the program output. This case has 
no mean flow, and the frequency is such that only one radial mode is cut-on for the m = 6 
spinning mode. The results will be compared below to a flow case. 


PROGRAM FDANC1 - TURBOFAN ANC SPECS 

CYLINDRICAL OR ANNULAR HARDWALL DUCT ROOTS FOUND AS 
ROOTS OF DERIVATIVE OF BESSEL FUNCTION 

INNER RAD = .000 OUTER RAD - 2.000 DUCT HGT = 2.000 ft 

DUCT MACH # = .00000 

FLOW TEMPERATURE = 70.00000 deg F 

SPEED OF SOUND = 1128.295 ft/sec 343.904 m/sec 

FREQUENCY = 935.0 Hz., ETA = 1.65737 

WAVE NUMBER, k = 17.08260 per m. , kH = 10.41355 

SPL = 110.000 dB 

AXIAL LENGTH OF ANC PORT = 4.500000 in. 

AIR DENSITY = 1.2300E+00 kg/rrT3 

RHO*CS = 423.0023 MKS Rayls 

SPINNING MODE ORDER, m = 6 

NUMBER OF RADIAL MODES = 6 



THIS DUCT TAKEN 

AS CYLINDRICAL 

WITH INNER 

RADIUS SET TO 

0 



JR 

HW ROOT 

PROP CNST 

(RE, IM) 

CUT-OFF FRQ 

HW COR 

AXL WVLNTH 

MOD ANG 


ALPHA* R2 

PER METER 

HZ 


ft 

deg 

1 

7.501266 

11.848900 

.000000 

673.515 

1.39 

. 5303 

46.08 

2 

11.734940 

.000000 

8.374456 

1053.643 

.89 

.0000 

90.00 

3 

15.268180 

-.000001 

18.316620 

1370.882 

.68 

.0000 

90.00 

4 

18.637440 

-.000001 

25.355620 

1673.397 

. 56 

.0000 

90.00 
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5 21.931720 - . 000001 31.663000 1969.180 .47 .0000 90.00 

6 25.183930 -.000002 37.614940 2261.186 .41 .0000 90.00 


PRESSURE AT SPL = 6.3246E+00 nt/nV-2 

ENERGY IN J=1 MODE = 2.7563E-02 watts 

DESIRED MODE AMPLITUDE = 1.7859E+01 nt/m^2 

INITIAL PISTON VELOCITY, INNER WALL = .000000 .000000 m/sec 

INITIAL PISTON VELOCITY, OUTER WALL = 1.000000 .000000 m/sec 

INNER VELOC/OUTER VELOC = .000000 .000000 


JR 

EV-ALPHA 


MODE C-VAL 

MODE COEFF ( RE, IM) 

MAG MODE AMPL 


PER METER 



nt/m'"' 

2 

nt/m'*' 2 

1 

12.305230 

.000000 

.000000 

1 . 7 8 58 88E+01 

-1 . 561 9O0E-07 

1 . 7 859E + 01 

2 

IS. 250220 

.000000 

.000000 

1 . 5354 35E-01 

1. 90991 6E-01 

2 . 4 506E-01 

3 

25.046230 

.000000 

.000000 

1 . 7 8 8268E-01 

2.22441 3E-01 

2 . 854 IE-01 

4 

30. 5*73230 

. 000000 

. 000000 

7 . 065713E-02 

8 . 788983E-02 

1.1277E-01 

5 

35 . 977220 

.000000 

. 000000 

8.507 2 93E-02 

1 . 058215E-01 

1. 3578E-01 

6 

41.312210 

.000000 

. 000000 

1 . 32567 8E-02 

1 . 64 9000E-02 

2 . 1158E-02 

VELOCITY AT INNER 

WALL * 

0 . OOOOOE+OO 

0. 00000E+00 m/s 



VELOCITY AT OUTER 

WALL = 

-5.77585E-02 

4 . 64 337E-C2 m/s 



DISPLACEMENT AMPLITUDE, INNER WALL = 

0. 00000E+00 in 




DISPLACEMENT AMPLITUDE, OUTER WALL = 


4.96668E-04 in. 


It should be noted that the values in the column marked “HW ROOT” are the 
roots of the Bessel function boundary condition, y m j, which are nondimensional, but that 
the value in the column marked “EV-ALPHA” are the oc,nj, in dimensions of 1 /meters. The 
propagation constants are also in units of 1/meters. The MODE C-VAL is the coefficient 
that multiplies the Bessel function of the second kind in the annular duct eigenfunction, 
and thus is zero for the cylindrical duct eigenfunction. The columns marked MODE 
COEFF (RE, IM) are the complex mode coefficients for the pressure solution, and are the 
desired results of the computation. 

Note that the desired mode amplitude to achieve 110 dB SPL in the first radial 
mode (the only cut-on mode) requires a pressure amplitude of 6.32 nt/m 2 at the outer wall 
and a mode coefficient amplitude of 17.86 nt/m 2 . The displacement amplitude that is 
required for the actuators to generate this mode coefficient is 5. 0x1 O' 4 in. 

The magnitudes of the mode coefficients for this case are plotted in Figure (3), and 
the pressure profile is shown in Figure (4). Note that the second highest mode, the third, 
is 36 dB down from the first mode, indicating excellent coupling to the cut-on mode for 
this case. The pressure profile is essentially the pure m = 6, j = 1 Bessel function of the 
first kind. 
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3.2.2 Computer Output for Case 2, Annular Duct, Mach 0.0, m = 6, Outer Wall 
Actuator Only 

This case is for the annular duct with no flow, with a radius ratio of 0.307. The 
actuator velocity at the inner wall in this case is set equal to zero. The output for the 
annular duct case is slightly different, and is shown below: 

PROGRAM FDANC1 - TURBOFAN ANC SPECS 

CYLINDRICAL OR ANNULAR HARDWALL DUCT ROOTS FOUND AS 
ROOTS OF DERIVATIVE OF BESSEL FUNCTION 

INNER RAD = .615 OUTER RAD = 2.000 DUCT HGT = 1.385 ft 

DUCT MACH # = .00000 

FLOW TEMPERATURE = 70.00000 deg F 

SPEED OF SOUND = 1128.295 ft/sec 343.904 m/sec 

FREQUENCY = 935.0 Hz., ETA = 1.14806 

WAVE NUMBER, k = 17.08260 per m. , kH = 7.21347 

SPL = 110.000 dB 

AXIAL LENGTH OF ANC PORT = 4.500000 in. 

AIR DENSITY = 1.2300E+00 kg/m A 3 

RHO*CS = 423.0023 MKS Rayls 

SPINNING MODE ORDER, m =* 6 

NUMBER OF RADIAL MODES = 6 


ANNULAR DUCT ANALYSIS USED 


J 

GREC 

GCYL 

GANN 





1 

. 0000 

7.5013 

7.5010 





2 

4 . 5353 

11.7349 

11.7178 





3 

9.0706 

15.2682 

15.1073 





4 

13.6059 

18.6374 

18.1882 





5 

18.1411 

21.9317 

21.5424 





6 

22 . 6764 

25.1839 

25.3502 





JR 

HW ROOT 

PROP CNST 

( RE , IM ) CUT 

’-OFF FRQ 

HW COR 

AXL WVLNTH 

MOD ANG 


ALPHA*R2 

PER METER 

H2 


ft 

deg 

1 

7 . 501023 

11.849320 

.000000 

673.493 

1.39 

. 5303 

46.08 

2 

11.717810 

. 000000 

8.813345 

1052.105 

.89 

. oooo 

90.00 

3 

15.107290 

-.000001 

17 . 954040 

1356.436 

. 69 

. 0000 

90.00 

4 

18 .188200 

-.000001 

24 .461990 

1633.061 

.57 

. 0000 

90.00 

5 

21 . 542440 

-.000001 

30.935500 

1934.228 

.48 

. 0000 

90.00 

6 

25.350170 

-.000002 

37 . 914260 

2276.112 

.41 

. 0000 

90.00 


PRESSURE AT SFL = 6.3246E+00 nt/nT2 

ENERGY IN J=1 MODE = 2.7591E-02 watts 

DESIRED MODE AMPLITUDE = 1.7859E-01 nt/nT2 

INITIAL PISTON VELOCITY, INNER WALL = .000000 .000000 m /sec 

INITIAL PISTON VELOCITY, OUTER WALL = 1.000000 . 000000 m/sec 


INNER 

VELOC/ OUTER 

VELOC = 

.000000 

. 000000 



JR 

EV-ALPHA 


MODE C-VAL 

MODE COEFF (RE, IM) 

MAG MODE AMPL 


PER METER 



nt/m"" 

2 

nt/m~2 

1 

12.304830 

-.000129 

.000000 

1 . 7 8588 8E + 01 

-1.703681E-07 

1 .7859E+01 

2 

19.222130 

-.014592 

.000000 

6.858777E-01 

8 . 531162E-01 

1 . 0946E+00 

3 

24.782310 

-.148402 

.000000 

6 . 013200E-02 

7.47 941 OE-02 

9 . 5 969E-02 

4 

29.836290 

-.451348 

.OOCOOO 

-7 . 54 5694E-02 

-9 . 38557 4E-02 

1 . 204 3E-01 

5 

35.338650 

-.392167 

.000000 

2 . 288 906E-02 

2. 847013E-02 

3. 6530E-02 

6 

41.564930 

.162785 

.000000 

-4 . 326065E-02 

-5 . 3808 96E-02 

6 . 904 3E-02 


VELOCITY AT INNER WALL = G.OOOCOE+OC 0 . OOOOOE+OO m/s 

VELOCITY AT OUTER WALL = -7.177S0E-02 S.77080E-02 m/s 

DISPLACEMENT AMPLITUDE, INNER WALL = O.OOOOOE+OO in 

DISPLACEMENT AMPLITUDE, OUTER WALL = 6.17243E-04 in. 
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Figure (5) shows the mode coefficient magnitude for this case, and Figure (6) 
shows the pressure profile. Note that the second radial mode is 24 dB down from the first 
mode in this case, indicating very good coupling efficiency. An actuator displacement 
amplitude of 6.2xl0' 4 inches is required to achieve the required SPL level. Generally, the 
case is very similar to the cylindrical duct case, as one might expect at this low radius 
ratio. 


Case 02 



Figure (5) Mode coefficients for Case 2 



Figure (6) Pressure profile for Case 2 
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3.2.3 Comparison of Cases 1 and 3 

Cases 1 and 3 are for the cylindrical duct, Case 1 with no flow and Case 3 with a 
flow Mach number of 0.1. The mean flow of Mach 0.1 has had very little effect on the 
mode coupling due to the actuator in this case, as can be seen from Figure (7), which 
compares the mode magnitudes for the cases with and without flow. The pressure profiles 
are nearly indistinguishable from those presented previously. 


Case 01 vs. Case 03 



Radial Mod* Humber 


Figure (7) Mode coefficients for Case 3 compared to coefficients for Case 1 . 


3.2.4 Comparison of Cases 2 and 4 

Cases 2 and 4 compare the annular duct with and without flow. Case 4 is the 
baseline case for the annular duct with the actuator active on the outer wall only, with 
Mach 0. 1 flow. As in the cylindrical duct case, the mean flow of Mach 0. 1 has had little 
effect. Figure (8) compares the modes generated by the actuator for Case 2 (no-flow) and 
Case 4 (Mach 0.1). The pressure profiles are nearly indistinguishable from those 
presented previously. 
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Case 02 vs. Case 04 


100.00 



1 2 3 4 5 8 

Radial Mode Number 


Figure (8) Comparison of mode coefficients for Cases 2 and 4. 


3.2.5 Comparison of Cases 4 and 5 

In Case 5, only the inner wall of actuators is active, otherwise the conditions are 
the same as for Case 4. Figure (9) compares the generated mode magnitudes for Cases 4 
and 5, and indicates that the actuators at the inner wall only tend to generate higher output 
in the higher order modes, an indication of reduced coupling efficiency. The second mode 
is 14.7 dB down in Case 5, compared to 24.7 dB down in Case 4. Figure (10) shows the 
slightly different pressure profile obtained from the inner wall actuator case. 
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problem. The objective is to choose the magnitude and phase of the inner actuators such 
that the combination of actuators gives the least amount of amplitude to the higher order 
modes, coupling most efficiently to the cut-on mode. There is no direct method of 
deriving this inner actuator velocity a priori (one would have to use the analysis developed 
in this study in the backward sense, rather than the forward sense of predicting what 
pressure modes couple to a given actuator motion). 

As a first estimate, the inner actuator velocity was set to a magnitude and phase 
determined by the ratio of the pressure first mode eigenfunction at the inner radius to the 
pressure at the outer surface, assuming the velocities might have roughly the same ratio. 
The inner to outer radius pressure ratio was determined to be 0.01 5432 + O.Oi. 

Figure (11) compares the mode coefficients for these two cases, and indicates that 
there is not much difference. The radius ratio in this case, however, is quite low, at 0.3 1, 
so that the velocity at the inner actuator is quite small. There may be some other value of 
inner wall actuator velocity that offers some improvement in coupling efficiency, but it 
would have to be sought by trial and error. 



Figure (11) Comparison of mode coefficients for Cases 4 and 6 


3.2.7 Comparison of Cases 7 and 8 

For Cases 7 and 8 the inner radius was increased to 1 .5 feet, giving an annular duct 
radius ratio of 0.75. In Case 7, only the outer wall actuator is activated. In Case 8, both 
inner and outer wall actuators are activated, and the ratio of inner wall actuator velocity to 
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outer wall is determined using the same method described in Case 6. Figure (12) shows 
the comparison of mode coefficients for Cases 7 and 8. 


Case 07 vs. Case 08 



Figure (12) Comparison of mode coefficients for Cases 7 and 8 


In this case, a marked improvement of coupling efficiency is noted. The second 
radial mode is 24.3 dB down in Case 7 with outer wall actuators only, but the second 
radial is 37 dB down for Case 8, for this particular choice of inner wall actuator velocity. 
Although this may not be the optimum coupling efficiency, it does indicate the potential 
for improved coupling to a single mode using both inner and outer actuator rings. 

The pressure profiles are almost identical in these cases, since the only difference is 
in the higher order mode content, which is small for both cases. It should be noted that 
the substantial 24 dB difference obtained with the outer wall actuator ring only may make 
the additional expense of an inner actuator ring unnecessary. It would certainly seem to 
be unwise for the low radius ratio case. This would have to be evaluated on a case-by- 
case basis. 


3.2.8 Discussion of Case 9 

Although the computer program is aimed primarily at the case in which only one 
radial mode is cut on, for which it calculates the “spillover” into cut-off modes, it will 
calculate the radial mode coupling independently of the number of cut-on radial modes. It 
is not obvious, however, how to choose the inner wall actuator velocity relative to the 
outer wall actuator velocity to best couple to two radial modes, which will have relative 
magnitudes and phases. It is not clear at this point whether the two coplanar actuator 
rings can be used to generate two independent radial modes. 
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As an example of a case with two cut on radial modes. Case 9 increases the 
frequency of Case 8 to 1350 Hz. The mode coefficients of the fan noise source are 
assumed to be the following values at the actuator source plane (assuming these modes 
are the anti-phase values): 

A, = 1.0 + 0.0 i 
A 2 = 0.5 -0.2i 

This is the source that the ANC actuator system must try to match. The difficulty is to 
determine the inner wall actuator velocity (relative to the outer wall), that will best match 
these two modes. 


In an actual ANC system with feedback control sensors, the inner and outer wall 
actuator velocities in the m = 6 mode would be adjusted (in somewhat of a trial-and-error 
process) by the control system to values that attempt to achieve cancellation. A key 
question is whether such a set of control parameters exist, because this will bear directly 
on the stability of the control system operation. 


In this case, a first guess at the inner wall velocity was made by computing the 
ratio of the pressure of the sum of the two fan source modes at the inner wall to the 
pressure of the sum of the two fan source modes at the outer wall. The following is the 
output run for this case: 


PROGRAM FDANC1 - TURBOFAN ANC SPECS 

CYLINDRICAL OR ANNULAR HARDWALL DUCT ROOTS FOUND AS 
ROOTS OF DERIVATIVE OF BESSEL FUNCTION 

INNER PAD = 1.500 OUTER RAD = 2.000 DUCT HGT = .500 ft 

DUCT MACH # = .10000 

FLOW TEMPERATURE = 70.00000 deg F 

SPEED OF SOUND = 1128.295 ft/sec 343.904 m /sec 

FREQUENCY = 1350.0 Hz., ETA = .59825 

WAVE NUMBER, k = 24.66471 per m. , kH = 3.75890 

SPL = 110. 000 dB 

AXIAL LENGTH OF ANC PORT = 4.500000 in. 

AIR DENSITY = 1.2300E+00 kc/m A 3 

RHO*CS = 423.0023 MKS Rayls 

SPINNING MODE ORDER, m = 6 

NUMBER OF RADIAL MODES = 6 


ANNULAR DUCT ANALYSIS USED 


J 

GREC 

GCYL 

GANN 

1 

. 0000 

7.5013 

6.8531 

2 

12.5664 

11.7349 

14.4285 

3 

25.1327 

15.2682 

26.0950 

4 

37.6991 

18.6374 

38.3452 

5 

50.2655 

21.9317 

50.7513 

6 

62.8319 

25.1839 

63.2209 


JR 

HW ROOT 
ALPHA* R2 

PROP CNST { RE, IM) 
PER METER 

CUT-OFF 

FRQ 

H2 

HW COR AXL WVLNTH 
ft 

MOD ANG 
deg 

1 

6.853070 

19.713180 .000000 

612. 

.231 

2.21 .3187 

26. 97 
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2 

14.428500 

4.913265 

.000000 

1288.996 

1.05 

1.2788 

72.71 

3 

26.095010 

-2.491387 

35.074630 

2331.245 

.58 

. 0000 

90.00 

4 

38.345170 

-2.491387 

58.102920 

3425.634 

.39 

.0000 

90.00 

5 

50.751250 

-2.491388 

79.877620 

4533. 954 

.30 

.0000 

90.00 

6 

63.220940 

-2.491389 

101.210100 

5647 . 956 

.24 

.0000 

90.00 


PRESSURE AT SPL - 6.3246E+00 nt/m^ 

ENERGY IN J=1 MODE = 4.9441E-02 watts 

DESIRED MODE AMPLITUDE = 1 . 9178E+01 nt/m~2 


INITIAL PISTON VELOCITY, INNER WALL = .280375 
INITIAL PISTON VELOCITY, OUTER WALL = 1.000000 
INNER VELOC/OUTER VELOC = .280375 .163253 


.163253 m/sec 
. 000000 in/sec 


JR EV- ALPHA 

PER METER 


MODE C-VAL 


MODE COEFF ( RE , IM ) 
nt/m^2 


MAG MODE AMPL 
nt/m A 2 


1 11.241910 

2 23.668810 

3 42.806770 

4 62.902170 

5 83.253360 

6 103.708900 


-.305718 

. 000000 

1. 917825E+01 

1 . 324236E-08 

1 . 9178E+01 

-. 960686 

. 000000 

-9. 313667E+00 

2 . 115661E+01 

2. 3116E+01 

1.222404 

. 000000 

4 . 968037E-02 

1 . 549457E-02 

5 . 204 IE-02 

. 355008 

. 000000 

6. 742401E-02 

9.3237 16E-03 

6. 8066E-02 

.063063 

.000000 

1 .237546E-02 

1 . SI 6386E-03 

1 . 2523E-02 

-.105775 

.000000 

2 . 008836E-02 

2 - 43067 3E-03 

2 . 0235E-02 


VELOCITY AT INNER WALL = 1.10319E-C2 -1 . 78015E-02 m/ s 

VELOCITY AT OUTER WALL = 1.77576E-C3 -6.45257E-02 m ,/s 

DISPLACEMENT AMPLITUDE, INNER WALL = 9.72090E-05 in 

DISPLACEMENT AMPLITUDE, OUTER WALL = 2.99620E-04 in. 


The output modes, normalized such that the first mode is unity, are given by 


A, = 1.0 + 0.0i 
A 2 = -0.486+ 1.1 03 i 

The second mode generated by the actuator rings does not closely match the desired fan 
source mode. A brief and cursory parametric study was made to attempt to find the inner 
wall actuator velocity that would match the desired mode coefficient, but this was 
unsuccessful in finding any improvement. 

Since this is just a single example case, it offers neither proof nor ability to 
generalize, but it may indicate that there will be problems attempting to match two cut-on 
modes with coplanar actuator rings. Further investigation of the possible advantages of 
inner and outer rings is warranted, and is recommended for future effort. 
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4. Conclusions and Recommendations 


4.1 Conclusions 

An analysis has been developed to predict the modes generated in an annular duct 
due to the coupling of flush-mounted ring actuators on the inner and outer walls of the 
duct. The analysis has been combined with a previous analysis for the coupling of modes 
to a cylindrical duct in a FORTRAN computer program FDANC to perform the 
computations. The method includes the effects of uniform mean flow in the duct. The 
program can be used for design or evaluation purposes for active noise control hardware 
for turbofan engines. 

Predictions for some sample cases modeled after the geometry of the NASA Lewis 
ANC Fan indicate very efficient coupling to the lowest order radial mode in both the inlet 
and exhaust ducts for the m = 6 spinning mode at frequencies where only a single radial 
mode is cut-on. Radial mode content in higher order cut-off modes at the source plane 
and the required actuator displacement amplitude to achieve 110 dB SPL levels in the 
desired mode were predicted. Equivalent cases with and without flow were examined for 
the cylindrical and annular geometry, and little difference was found for a duct flow Mach 
number of 0. 1 . 

One case using both inner and outer actuator rings in a high-radius-ratio annular 
duct indicated that, when only one radial mode is cut-on, it might be possible to phase 
both rings to give a more efficient match to the desired mode than using only the outer 
ring. Whether the increased matching efficiency would justify the added hardware 
expense, however, is questionable. 

One case was considered at a higher frequency in which two radial modes were 
cut-on, in a preliminary attempt to generate anti-phase matches to both modes using 
coplanar inner and outer actuator rings. In brief search, no combinations of inner and 
outer actuator velocity that matched the given modes could be found. Although the 
search was far from extensive, this may indicate that two coplanar rings may have a 
difficult time matching two arbitrary duct modes. 


4.2 Recommendations for Further Effort 

Further checkout and de-bugging of the prediction program is required to achieve 
a high level of confidence in its validity. This would include some parametric studies in 
which frequency, radius ratio, and flow Mach number are varied. It is recommended that 
the analysis be checked against test data obtained from ANC testing in the NASA Lewis 
4-foot Fan Rig. 

The actuator ring coupling program can be adapted as a subroutine to the 
cylindrical duct modal analysis and the exhaust duct modal analysis. This will allow the 
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fan source to be defined in terms of characteristic modes at the fan source plane and 
predict the propagation to the arbitrarily-located ANC source plane. The actuator 
velocities can then be determined to generate the anti-phase mode. The resulting 
combined fan source/ ANC pressure can then be calculated at any desired wall sensor 
position. The actuator velocities can be determined manually or using a simulation of a 
control system feedback loop. 

It would be useful to determine effects of higher Mach numbers on mode coupling, 
since this preliminary study went only to Mach 0.1 (representative of the flow in the 
NASA ANC fan). For the case with one radial mode cut-on, it would be useful to develop 
an optimization scheme that determines the inner wall actuator velocity that gives best 
coupling to the radial mode, and determine the implications of this increase in efficiency in 
terms of required power input, etc. 

Finally, the possibility and limitations of coupling to two cut-on radial modes 
using coplanar inner and outer actuator rings should be examined. A plot of the mapping 
of the complex inner wall actuator velocity to the complex mode coefficient of the second 
radial mode could be made for specific cases The implications of these results to design 
of feedback control systems should be explored. 
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5. Appendices 

5.1 Appendix A - Numerical Computation of Bessel Function Integrals 

The solution of the non-homogeneous differential equation that arises from the 
problem of the wall velocity boundary condition in the cylindrical duct requires the 
computation of four integrals involving Bessel functions of the first and second kind in the 
integrands. The integrals requiring evaluation are: 


= J(Pmn(amn r ) dr 

r i 


(96) 


Hr = j (Pmn( a nmr) rdr ( 97 ) 

f l 


Hr = J <Pnrn( a mnr) r2dr ' ( 98 ) 

r i 


and 


where 


Hr = J(Pmn(anmr) r3dr 


( Pmn( a mn r ) ^m( a mn r ) ^mn^m( a nin r ) 


(99) 


( 100 ) 


It may be possible to reduce the order of the Bessel function in the integration by 
using identities such as 


J x n (p mn (x)dx = -x n cp m .i, n ( x ) + ( n + m-i)j x n ' 1 (p m .] n (x)dx 


( 101 ) 


and, for certain m-values, this will reduce to a closed form in Bessel functions. Generally, 
however, one is left with the integral of a Bessel function that must be computed 
numerically, since there is no general integration formula for the integral 

Jcp mn (x)cix 
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where m is an arbitrary integer. Since one eventually ends up with an integral to be 
evaluated numerically, it seems more straightforward to simply integrate the original forms 
of the integrals by numerical integration. 

A series of numerical integration subroutines have been adapted to evaluate these 
integrals. The subroutines are written in FORTRAN, suitable for use on the EBM-PC, and 
can be easily adapted to almost any other platform. The integrals are solved numerically 
by the Romberg integration method based on subroutines provided in Numerical Recipes'. 
The subroutine CBESEL, developed at GEAE, is used to compute the Bessel functions. 
The subroutines have been checked out against a MathCAD computation. 


5.2 Appendix B - Contour Integrations for Inversion of Fourier Integrals 

The Fourier Transform Inversion Integral required to find the z-dependent modal 
coefficients of the solution to the wall velocity boundary condition equation can be split into six 
separate integrals, which, in general form, are 


Gj(z) = 
forj = 1, 2, or 3, or 
Gj(z) = 

for j = 4, 5, or 6, where the Q™ are constant coefficients that depend on m and n. %{ yl ) is the 
Fourier Transform of the axial velocity function for the ANC noise source. 

Following Butkov 5 , these integrals can be evaluated as contour integrals in the complex 
plane using the theory of residues. The path of integration must be prescribed, taking into 

account poles in the complex plane. Assuming that the velocity function transforms K(k) or 

H(k) have no zeros or poles in the region surrounded by the contour, the poles of the 
integrands of both integrals are at the complex points given by the roots of the denominator in 
the integrand: 


Q(j) 
V mn 


I 


K 


j-1 


%{ K)( 


-1KZ 


2% ^ (k - aL) + 2kMK - (l - M 2 )k^ 


dK 


( 102 ) 


q(j) 

Vmn 

271 


J 


HjOOe 


-IKZ 


(k - a^) + 2kMx - (l - M 2 )k : 


dK 


(103) 


Jj) -kM±^k 2 -(l-M 2 )aL 
K 1 - M 2 


( 104 ) 


Figure (13) shows the contour, which is chosen in a particular manner that results only 
in waves traveling in the positive z-direction. The contour consists of the negative and 


37 



positive real axis and a semi-circle of infinite radius that is closed downward. The integration 
direction is counter-clockwise. Since the semi-circular part of the integral at infinite radius 
doesn't contribute in the lower half plane {e 1KZ -» 0), the integral, by the residue theorem, is 2%i 
time the sum of the residues enclosed by the contour. 



Figure (13) Shape of the contour for contour integration. 
Thus 


Gj(z) 


2rciQ ( njn Res [. _L X(K)e iB 

1-M 2 _ ( K - k^)(k - k^)_ k = ic (2) 


( 105 ) 


and, since the residue at k = k ( 2) is 


Lim 

k-»k( 2 > 


1 ff(K) e - iKZ ~ 

K - 



( 106 ) 


the first integral becomes 


G](z) 





( 107 ) 


Following the same reasoning and using the same contour for the second integral, it 
becomes 
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( 108 ) 


Gj(z) = - 27iiQ^Res 


1 K (j ~ 1) ^(K)e~ iKZ 


2k (k - k (1) )(k - k (2) ) 


IK = 


K 


( 2 ) 


which, by the residue calculus, reduces to 


Gj(z) = i(ic 


+ 

mn 


V^qO) 

f Vmn 


%< 2 >) 


JK^Z 


^k 2 -(l-M 2 )a 


2 

mn 


( 109 ) 


Making the correct substitutions for k (2) and Q™, one can use these formulae to evaluate all 
required inverse transforms. The formulae apply to the cases with EU, H 5 , and H« if these 
transforms are sufficiently well-behaved, otherwise they must be evaluated as special cases. 
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